Cold Atoms and Molecules in Self- Assembled Dipolar Lattices 
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We study the realization of lattice models, where cold atoms and molecules move as extra particles 
in a dipolar crystal of trapped polar molecules. The crystal is a self-assembled floating mesoscopic 
lattice structure with quantum dynamics given by phonons. We show that within an experimentally 
accessible parameter regime extended Hubbard models with tunable long-range phonon-mediated 
interactions describe the effective dynamics of dressed particles. 



Trapped atomic and molecular quantum gases allow 
the realization of quantum lattice models of strongly in- 
teracting bosonic and fermionic particles. For example, 
the dynamics of atoms in optical lattices is well described 
by a Hubbard model, where the tunability of the Hub- 
bard parameters via external fields combined with atomic 
physics techniques of preparation and measurement pro- 
vides a quantum simulator of strongly correlated con- 
densed matter models In this letter we propose and 
study an alternative scenario of realizing lattice models, 
where a dipolar crystal of trapped polar molecules pro- 
vides a self-assembled floating lattice structure for ex- 
tra particles, which are atoms or molecules of a second 
species (Fig. 1). By confining polar molecules to an ef- 
fective 2D (Fig. la) or ID (Fig. lb, c) geometry by strong 
transverse trapping, dipolar crystals can form as a result 
of the balance of strong repulsion between the dipoles 
aligned by an external electric field, and an in-plane trap- 
ping potential 0, 0] • Particles moving in this lattice un- 
der conditions of elastic scattering see a periodic poten- 
tial, and thus form a lattice gas. 

The distinguishing features of this realization of lattice 
models are: (i) Dipolar molecular crystals constitute an 
array of microtraps with its own quantum dynamics rep- 
resented by phonons (lattice vibrations) , while the lattice 
spacings are tunable with external control fields, ranging 
from a /xm down to the hundred nm regime, i.e. poten- 
tially smaller than for optical lattices, (ii) The motion 
of the extra particles is governed by an interplay of Hub- 
bard (correlation) dynamics in the lattice and coupling 
to phonons. The tunability of the lattice allows to access 
a wide range of Hubbard parameters and phonon cou- 
plings. Compared with optical lattices, for example, a 
small scale lattice yields significantly enhanced hopping 
amplitudes, which set the relevant energy scale for our 
Hubbard model (e.g. for exchange interactions), and thus 
also the temperature requirements for realizing strongly 
correlated quantum phases. While the molecular setups 
we propose are reminiscent of (and may be relevant to) 
solid state systems with strong phonon couplings, as e.g. 
polaronic and/or superconducting materials, they real- 
ize the unusual parameter regime where the mass of the 
crystal's and of extra particles can be comparable. 

A homogeneous lattice of polar molecules underlying 
the configurations of Fig. 1 relies on the strong repul- 
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Figure 1: A dipolar crystal of polar molecules in 2D (a) and 
ID (b,c) provides a periodic lattice V cp for extra atoms or 
molecules giving rise to a lattice model with hopping J and 
long-range interactions Vi,j (see text), (a) In 2D a triangu- 
lar lattice is formed by polar molecules with dipole moment 
d c perpendicular to the plane. A second molecular species 
with dipole moment d p <C d c moves in the honeycomb lattice 
Vcp (darker shading corresponds to deeper potentials), (b) A 
ID dipolar crystal with lattice spacing a provides a periodic 
potential for a second molecular species moving in a paral- 
lel tube at distance b (Setup 1). (c) ID setup with atoms 
scattering from the dipolar lattice (Setup 2). 



sive dipole-dipole interactions 14(11) = d'c/R 3 with R 
the distance between the molecules, and d c the dipole 
moment induced by a transverse electric field E^c- A 
requirement for the existence of a crystal is that the 
ratio of the potential energy to the kinetic energy of 
small oscillations around the equilibrium position, r<2 = 
d 2 m c /h 2 a, with a the lattice spacing and m c the mass, 
is larger than a critical value r c , where r c = 18 ± 4 
and r c ~ 1 for bosons at zero temperature in 2D and 
ID, respectively 0|. Thus a dipolar crystal will form 
for a < a max = d 2 m c /h 2 r c . In addition, we have 
flmin < a with a m i n = {\2d 2 /m c Lo' 1 _) 1 / b , which reflects 
the requirement of strong transverse trapping with a 
harmonic oscillator frequency u>± to prevent collapse 
due to attractive forces between aligned dipoles. For 
RbCs (SrO) molecules with permanent dipole moment 
d c = 1.27D (d c = 8.9D) confined by an optical lattice 
with uj±/2it ~ 150kHz, a m - m ~ lOOnm (200nm), while 
flmax can be several /mi. Excitations of the crystal are 
acoustic phonons with Hamiltonian H c = Y^ q ^ UJ q a l a q^ 
where a q destroys a phonon of quasimomentum q in 
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the mode A. In ID, Hu q = (2/tt 2 ) [12r d / 9 ] 1/2 £ RjC 
with Debye frequency Huju = Hcu^/a ~ l-4y/rXE/R )C , 
/« = Sj>o 4sin(gaj/2) 2 /j 5 , and lattice recoil frequency 
En c = h 2 ir 2 /2m c a 2 (typically a few to tens of kHz). The 
classical melting of the crystal into a normal phase oc- 
curs at ksTc — 0.018rdi?R :C and 0.2r £ j£ , R iC in 2D and 
ID, respectively 0,0]. 

An extra particle confined to the 2D crystal plane (Fig. 
la) or a ID tube (Fig. lb,c) will scatter from the periodic 
lattice potential J2j Kp(Rj' — r) with r and Rj the coordi- 
nates of the particle and crystal molecule j, respectively. 
We write = R°+Uj with R° the equilibrium positions 
and Uj small displacements, assuming that the particles 
do not significantly perturb the lattice. For particles be- 
ing molecules, this potential is given by the repulsive 
dipole-dipole interaction V cp (Rj — r) = d p d c / \ R/ — r | 3 
with dp <§C d c the induced dipole moment, and in the case 
of atoms we assume that the interaction can be modeled 
by a short range pseudopotential proportional to an elas- 
tic scattering length a cp . In addition, extra molecules and 
atoms will interact according to dipolar, or short range 
interactions, respectively. 

We consider a situation where the dynamics of the ex- 
tra particles in the lattice can be described by a sin- 
gle band Hubbard Hamiltonian coupled to the acoustic 
phonons of the lattice 0| 

H = -J J2 4 + ^^^ 

<i,j> i,j 

+ ^■ ! <' <vH ''Yj K + oL,) + H c . (1) 

1,3 

The first line describes the nearest neighbor hopping of 
the extra particles with hopping amplitudes J, and in- 
teractions V. We denote by a (c\) destruction (creation) 
operators of the particles. The first term in the sec- 
ond line is the phonon coupling obtained in lowest or- 
der in the displacement Uj = i Y q {h/2m c Nw q ) 1 l 2 £, q {a q + 

a%)e tq R ? with M q = F q q • ^(^/2A^m c u; g ) 1 / 2 ^ q , where 
£ g and TV are the phonon polarization and the num- 
ber of lattice molecules, respectively, Vq is the Fourier 
transform of the particle-crystal interaction V cp , and 
/3 q = J dr|wo(r)| 2 e lqr , with wo(r) the Wannier function 
of the lowest Bloch band [f| . The validity of the single 
band Hubbard model requires J, V < A, and tempera- 
tures ksT < A with A the separation to the first excited 
Bloch band. We note that the Hubbard parameters are of 
the order of magnitude of the recoil energy, J, V ~ -Er,c, 
and thus they are (much) smaller than the Debye fre- 
quency hw D ~ E R ,cy/rl, for r d > 1 @|. 

Below we will present detailed results for the examples 
of Figs. lb,c. The separation of time scales J, V-^i fou)D, 
combined with the fact that the coupling to phonons is 
dominated by high frequencies hw > J, V (see the dis- 
cussion of M q below) is reminiscent of polarons as parti- 
cles dressed by (optical) phonons, where the dynamics is 
given by coherent and incoherent hopping on a lattice 0, 



0] • This physical picture is brought out in a master equa- 
tion treatment within a strong coupling perturbation the- 
ory. The starting point is a Lang-Firsov transformation 
of the Hamiltonian H — > SHS^ with a density-dependent 
displacement S = exp[—^ q j(M q /h^ q )e lclIi 3 CjCj(a q — 

This eliminates the phonon coupling in the sec- 
ond line of Eq. (UJ in favor of a transformed kinetic 
energy term — JY2<% j> c i c j-^i-^j> where the operators 
Xj = exp[^ g M g e iqR 3 [a q — aL q )/tko q ] can be inter- 
preted as a lattice recoil of the dressed particles in 
a hopping process. In addition, the bare interactions 
are renormalized according to = + with 

V}p = -2£ 9 cos(q(R° - R°))M 2 //^ ? , that is, the 
phonon couplings induce and modify off-site interactions. 
The onsite interaction is given by Vjj = Vjj — 2E p 
with E p = ~^2 q M 2 /hii) q the polaron self-energy or po- 
laron shift. For J — the new Hamiltonian is di- 
agonal and describes interacting polarons and indepen- 
dent phonons. The latter are vibrations of the lattice 
molecules around new equilibrium positions with un- 
changed frequencies. For the models of Fig. 1 b,c stud- 
ied below, consistency with the assumption of small per- 
turbation of the lattice by the extra particles requires 
Au/a = £ 9 (8/i 2 /m c N) 1 / 2 M q sm{qa/2) 2 /{hw q f/ 2 a < 1, 

which is achieved for {E p E-r. c ) x / 2 <C hujo, and is gener- 
ally satisfied. 

A Born-Markov approximation with the transformed 
kinetic energy as perturbation, and the phonons a finite 
temperature heatbath with J, V -C fnujj (see above), pro- 
vides us with the master equation for the reduced density 
operator of the dressed particles p(t) in Lindblad form, 

m = £ &ffb jS b w ,p(t)} (2) 

j,l,SS> 

+ -^-( 2b jSP(t)biS' - b jS bis>p(t) ~ p^bjsbig'), 

j,l,S,S> 

with bjs = Cj +g Cj, and where the effective system Hamil- 
tonian becomes an extended Hubbard model, 

H = -jJ2ctc j + ^yiAc}o j c i . (3) 
<*>i> i,3 

Coherent hopping of the dressed particles is described 
by J = J((XjXj)) = Jexp(-S T ), where S T = 
Y^q {M q /hijj q ) 2 [1 — cos(qa)] (2n q (T) + 1)] characterizes the 
strength of the particle-phonon interactions, and n q (T) is 
the thermal occupation at temperature T [5|]. The terms 

involving A^'f are second order corrections, which are 

small relative to H both in the "weak" St=o = So < 1 
and "strong" So ^> 1 coupling regimes, provided J <C 
hwo and J -C E p , respectively [§]. The dissipative 

term in Lindblad form in the second line of Eq. ([2]) cor- 
responds to thermally activated incoherent hopping with 
rates Tj'f , which are small compared to J for St -C 1 
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and S T > 1, provided JV 2 =0 k^mw D ) 4 ^\ < 1 
and ksT/Ep <C 1, respectively 0, @|, and in particu- 
lar they are negligible for the energies of interest fcgT <C 
min(A, E p , k B T c ). 

In the parameter regime of interest the dynamics of the 
dressed particles is described by the extended Hubbard 
Hamiltonian H. In the following we calculate the effec- 
tive Hubbard parameters from the microscopic model for 
the ID setups described in Figs, lb and c [1(J |. 




Figure 2: Setup 1 (Fig. lb): Hubbard parameters for d p /d c = 
0.1 and m c = m p . (a) Lattice depth Vb in units of i?R,c vs. b/a 
for Td = 50 and 500. Thick continuous lines: tight-binding 
regime 4 J < A. (b) Reduction factor So (dashed dotted lines) 
and polaron shift E p /J (solid lines), for 4 J < A. (c) Contin- 
uous lines: phonon- mediated interactions Vjj +1 . Horizontal 
(dashed) lines: Vj,j+i. (d) Contour plot of Vjj+i/2J (solid 
lines) as a function of b/a and ra. A single-band Hubbard 
model is valid left of the dashed line (4 J, Vij < A), and right 
of the dot-dashed line (E p < A). 



In Setup 1 (Fig. lb) molecules of a second species 
are trapped in a tube at a distance b from the crys- 
tal tube under ID trapping conditions. For crystal 
molecules fixed at the equilibrium positions with lat- 
tice spacing a, the extra particles feel a periodic po- 
tential V cp (x) = d c d p J^j [b 2 + ( x — i a ) 2 \ 3 ^ 2 , which de- 
termines the bandstructure. The potential is sinusoidal 
for b/a > 1/4, while for b/a < 1/4 it has a comb- 
like structure, since the particles resolve the individ- 
ual molecules forming the crystal. The lattice depth 
Vb = V cp (a/2) — V cp (0) is shown in Fig. 2a as a function 
of b/a, where the thick solid lines indicate the parame- 
ter regime 4 J < A. The strong dipole-dipole repulsion 
between the extra particles acts as an effective hard-core 
constraint [ll|. We find that for 4 J < A and d p <C d c the 
bare off-site interactions satisfy Vij ~ d 2 /(a\i — j\) 3 < A. 



The particle-phonon coupling is 

with /Ci the modified Bessel function of the second kind, 
and M q ~ Jq for q — > 0. For b/a < 1, which is the 
regime of interest (compare Fig. 2) , M q is peaked at large 
q ~ n/a, so that the main contribution to the integrals 
in the definition of St and E p is indeed dominated by 
large frequencies hw q > J. A plot of So as a function 
of b/a is shown in Fig. 2a. We find the scaling So cx 
y/r2(d p / d c ) 2 , and within the regime of validity of the 
single band approximation, Sq can be tuned from So <C 
1 (J ~ J) to So 3> 1 (J <C J) corresponding to the 
large and small polaron limit, respectively. The polaron 
shift E p , as shown in Fig. 2b, generally exceeds the bare 
hopping rate J, and in particular, E p J for So > 1, 
Together with the condition hwu ~3> J this ensures that 
the second-order shifts in the master equation, A^ , are 
indeed negligible in the parameter regime of interest. 

The phonon-mediated interactions vf/j show oscilla- 
tions, which for b/a < 1/4 decay slowly as ~ l/\i — j\ 2 . 
These interactions are thus long-ranged, and, depend- 
ing on their sign, they can enhance or reduce the direct 
dipole-dipole repulsion of the extra particles. The term 
vjj +1 is shown in Fig. 2c to alternate between attrac- 
tive and repulsive as a function of b/a. The effective Hub- 
bard parameters Vjj+i and J are summarized in Fig.HJl, 
which is a contour plot of Vj t j+i/2J as a function of 
and b/a. The ratio V/, J+ i/2J increases by decreasing 
b/a or increasing rj, and can be much larger than one. 
Equation (J3J) is valid in the region left of the dashed line, 
where 4J, V < A, and right of the dashed-dotted line, 
where E p < A. For E p > A a multi-band approach is 
required. 

In Setup 2 (Fig. lc) neutral atoms are trapped in the 
same tube as the crystal molecules. For fixed molecules 
providing the lattice structure each atom feels the ID 
potential V cp (x) = g cp S(x — ja), which determines 
the bandstructure (Kronig-Penney model). Here g cp = 
2h 2 a cp /fi p a 2 L p for a 3D scattering length smaller than the 

transverse confinement, a rn <C a -L,p = (fi/mpWi.p) 1 ' 2 , 
with n the reduced mass [13]. For g cp /aE^ p ^> 1 the 
width of the lowest band is 4J ~ (pcp/ai^p) -1 , while 
the gap A tends to A ~ 3-Er, iP , with -Eh. iP = h 2 Tr 2 /2m p a 2 . 

In the following we are interested in bosonic atoms in- 
teracting with each other via a contact potential with 
coupling strenght g pp determined by their 3D scattering 
length cipp, which is tunable independent of a cp . The 
bare Hubbard interactions are dominated by onsite in- 
teractions, Vij ~ SijU, where U = 2>g pp /2a for large 
g cp . For the choice of parameters of Fig. 3 a descrip- 
tion in terms of a single-band Hubbard model is valid for 
g cp /aE R . p > 0.5 so that J, U < A. 
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Figure 3: Setup 2 (Fig. QJ): Hubbard parameters for 
m p /m c = 0.814, a pp = 2nm and uj± tC /2ir = 200kHz. (a) Re- 
duction factor So (dashed lines) and polaron shift E p /J (solid 
lines) as a function of the particle-phonon coupling strength 
g/aEn.p- The tight binding region, 4J, Vj,j < A, is right of 
the shaded area, (b) Effective onsite Vj,j/J (solid lines) and 
nearest-neighbor Vj,j+i/J (dashed lines) parameters. 



The particle-phonon coupling is 



2h 



M q = 

a V Nm c LU a 
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which is peaked at large q ~ ir/a, and M q ~ ^fq 
for q — > 0. We find the reduction factor So ~ 

0. 92(g cp /ai?R,p) 2 (m c /TOp) 2 /ry 2 and the polaron shift 
£" p — 0.94£;R ; p(5 C p/aSR iP ) 2 (TO c /mp)/r d , which decrease 
with increasing (see Fig. 3a). For the parameters 
of Fig. 3 we have So < 1 over a wide range of r^, 
so that J ~ J. The strong coupling regime So 3> 1, 
with Ep/J » 1, can be reached by decreasing rd or 
increasing g C p/aE Rp . We note that g C p/aE Rp is re- 
stricted by the condition of a stable crystal Ait <C a, 

1. e. g cp /aE R , p < 3r d . 

While the bare atom-atom interaction provides only 
an onsite shift, the phonon coupling induces long-range 
interactions which decay as ~ l/\i — j\ 2 . The ef- 
fective interactions Vj t j, Vj tt j+i and J are summarized in 
Fig. 3b, where the nearest-neighbor term Vjj+i ~ 0.16-E p 



is shown to be repulsive for all values of g cp /aE RiP , while 
the onsite interaction Vjj = U — 2E p turns from posi- 
tive to negative, which for bosons indicates an instabil- 
ity towards collapse. That is, for stability we require 
g pp /aE R>p > lA{g cp /aE K , p ) 2 (m c /m p ) 2 /rd. Fig. 3b 
shows that a regime of strong interactions Vjj,Vjj + i » 
J can be reached for r d ^> 1. 

One feature of our extended Hubbard model is the ap- 
pearance and tunability of strong off-site interactions, 
a necessary ingredient for a variety of new quantum 
phases (15 , [T3 |. which is difficult to realize in a stan- 
dard atomic setup As an example, at half fill- 
ing the particles in Setup 1 undergo a transition from 
a (Luttinger) liquid (V^i+i < 2J) to a charge-density- 
wave (CDW) > 2J), which can be observed, e.g. 
for r d = 100 at b/a « 0.5 (see Fig. 2d) p3|. Sim- 
ilarly, the groundstate of Setup 2 at half filling is a 
CDW for V jtj ,V jtj+1 > 4 J, pjj]. For larger filling and 
strong interactions 2Vjj+i > Vjj > 4J, the system un- 
dergoes a second order transition to a supersolid phase, 
where diagonal and off-diagonal orders coexist. These 
strong interactions are here realized, e.g., for rd = 250 at 
g cp /aE RtP ?=s 10, where Vjj+i « Vj.j w 7. 6 J (see Fig. 3b). 

In conclusion, we have studied a scenario where cold 
atoms or molecules move in the periodic potential pro- 
vided by a dipolar molecular crystal, with quantum dy- 
namics given by phonons. This is in contrast to familiar 
traps in atomic physics, where backaction is negligible. 
Strong phonon-mediated off-site interactions and particle 
localization open new perspectives for studying the inter- 
play between strong correlations and phonon dynamics in 
a tunable setup. 
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